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Abstract 

We analyze spin-dependent carrier dynamics due to incoherent electron-phonon scattering, which 
is commonly referred to as Elliott-Yafet (EY) spin-relaxation mechanism. Eor this mechanism one 
usually distinguishes two contributions: (1) from the electrostatic interaction together with spin¬ 
mixing in the wave functions, which is often called the Elliott contribution, and (2) the phonon- 
modulated spin-orbit interaction, which is often called the Yafet or Overhauser contribution. By 
computing the reduced electronic density matrix, we improve Yafet’s original calculation, which 
is not valid for pronounced spin mixing as it equates the pseudo-spin polarization with the spin 
polarization. The important novel quantity in our calculation is a torque operator that determines 
the spin dynamics. The contribution (1) to this torque vanishes exactly. From this general result, 
we derive a modified expression for the Elliott-Yafet spin relaxation time. 

TAGS numbers: 75.78.Jp, 72.25.Rb, 76.30.Pk 
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I. INTRODUCTION 


Spin relaxation plays a role in spin-dependent dynamics both on long and short timescales. 
While its contribution to (precessional) magnetization damping is usually ascribed to spin- 
lattice relaxation and treated phenomenologically, the analysis of ultrafast demagnetization 
dynamics has often been based on the microscopic concept of spin-flip processes due to 
electron-phonon interactions as developed for semiconductors in the 1950s.” Overhauser- 
was the hrst to identify the modulation of the spin-orbit interaction by lattice vibrations as 
the microscopic mechanism for spin relaxation due to incoherent electron-phonon scattering 
processes. Elliott^ argued shortly thereafter that there is an additional contribution to the 
spin relaxation due to the momentum-dependent spin mixing in the wave functions and 
that, consequently, even spin-diagonal incoherent scattering processes due to spin-diagonal 
interactions can contribute to spin relaxation. The subject was taken up again by Gerasi¬ 
menko and Andreev,” and Yafet.- The latter calculated the spin-flip matrix element due 
to electron-phonon interaction (as a function of electronic momentum transfer q) including 
both the Overhauser and the Elliott contributions and showed that the hrst few orders in q 
vanish due to a cancellation of Overhauser and Elliott contributions. Nowadays, this com¬ 
bination of Overhauser and Elliott contributions is usually called Elliott-Yafet mechanism 
because Yafet derived a relatively simple result for the close-to-equilibrium spin relaxation 
time, which is suitable for evaluation from ab-initio input and can be used to justify an 
approximate relation between the spin relaxation time and the momentum relaxation time. 
At present, both simplihed and ab-initio based expressions due to the EY mechanism are 
widely used, in particular for magnetization dynamics in metals.-'— 

In this paper, we present a new analysis of the spin-relaxation problem due to incoherent 
electron-phonon scattering, as it was originally considered by Overhauser, Elliott, and Yafet. 
We derive the dynamical equation for the change of the reduced electronic spin density 
matrix by expressing the spin dynamics in terms of phonon assisted density matrices. This 
approach achieves a correct description of the dynamics of the spin vector, as opposed to 
the Yafet treatment, which gets spins and pseudospins mixed up, and thus cannot correctly 
account for the “amount of spin-flip” in each scattering transition. If one correctly describes 
the spin vector, it becomes obvious that the important quantity for spin dynamics is a torque 
matrix element, which is not present in the conventional derivation. - Remarkably, there is no 
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contribution to this torque matrix element from spin-diagonal scattering mechanisms. Put 
differently, spin-diagonal electron-phonon scattering and spin-orbit coupling alone, which is 
usually referred to as the Elliott spin-relaxation mechanism, yields no spin dynamics. Based 
on this observation, we derive a modihed result for the EY spin relaxation time. 

This paper is organized as follows. In Sec. |TT]we briefly review the conventional Elliott- 
Yafet treatment, some more recent contributions, and how the Elliott-Yafet spin relaxation 
mechanism has been applied in theoretical models for the demagnetization of ferromagnets. 
In Sec. mil we set up the electron-phonon interaction hamiltonian, discuss long and short- 
range contributions, and derive the equations of motion for the spin density matrix and the 
spin expectation value. Section HVl is devoted to the derivation of a spin relaxation time for 
the special case of Kramers degenerate bands. An important ingredient for this derivation is 
the form of the quasi-equilibrium spin-density matrix in the presence of spin-orbit coupling, 
which is discussed in some detail. The conclusions are presented in Sec. |Vl and appendix lAI 
contains a short demonstration concerning the form of the spinor wavefunctions on Kramers 
degenerate bands. 


II. ELLIOTT-YAFET APPROACH 


The Elliott-Yafet approach has been reviewed often.->ii>i^ Only for comparison with our 
calculations, we repeat here some of the main points of Yafet’s derivation in Ref. [s using 
his notation. The objective of Yafet is to calculate the rate of a spin-flip transition for two 
Kramers degenerate bands including spin-orbit coupling. The Kramers degeneracy implies 
that for a band index b we have two wave functions 


V'b^k(x) = 06t,k(x)| t) + XH,k(x)| i) 

(1) 


(2) 


with the same energy Ebk, where K is the time-reversal operator. Focusing on one band and 
dropping the corresponding b index he then calculates the Golden-Rule transition probability 
lK|^k,jj.k' for a spin-flip transition —)■ due to the electron-phonon interaction. More 

precisely, Yafet calculates the dynamics of the spin polarization, which is defined as 


T^spin ^^[r^il"k rijj,k;]. 


( 3 ) 
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where riAk denote the carrier distributions with momentum k and pseudospins A =i]', 'll- It 
is clear that -Dgpin cannot be a good approximation to the spin polarization for pronounced 
spin mixing. Then (including a factor of 2 as in Yafet’s derivation) one obtains 


d 

dt 


Apin - - IA,ir) - 


( 4 ) 


Thus the change of spin polarization is determined essentially by the number of transitions 
(per unit time), which is obtained by adding the in-scattering and out-scattering Golden- 
Rule probabilities 


^ |lAk,ltk' %k[l — %k'] — IAk,.l^k' %k[l — %k']|, 

k.k' 


( 5 ) 


where 


27r, 


rW 


IRfrk^k' ^ l^^k,^k'| 2 MNu, 


5(A ~ Ek' + EjJq)Nq -|- 5(A ~ A' ~ EjJ-q){N-q -\- 1) (6) 


and q = |k' — k|. The other symbols have an obvious meaning and are defined below. To 
obtain a relaxation time valid for a small spin polarization, the distribution functions riAk are 
assumed to be of the quasi-equilibrium form /(e^k —A^a) with qb\ the (pseudo)spin dependent 
chemical potentials, and the distributions are expanded for small This treatment has 
recently been extended to ferromagnets.— Finally, Yafet shows from a symmetry argument 
that there is a cancellation between the two contributions to the electron-phonon (e-pn) 
coupling matrix element in the long-wavelength limit. In this paper, we refer to these two 
contributions as the Elliott and Overhauser contributions, which we explain in detail below. 

In our view this method has three problems. 


1) The spin polarization iAgpin, defined in ([3]), is computed as the difference of pseudospin 
occupation numbers as if the electrons were in pure spin states for all k values. Here, 
the Yafet derivation ignores spin mixing, which leads to a k dependent expectation 
value of Sz- Importantly, the modulus of this spin expectation value may be signifi¬ 
cantly smaller than h/2. 

2) The calculation of spin dynamics is based on transition probabilities between pseudo¬ 
spin states or, equivalently, pseudo-spin occupation numbers n.(^k and nj^k- Such 
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a treatment neglects coherences between the pseudospin states, which are the off- 
diagonal elements of the spin density matrix 


P(k) 




(7) 


3) Yafet proves the cancellation between the Elliott and Overhauser contributions to 
the spin-flip matrix element for the short range part of the electron-phonon 

coupling matrix element. As shown by Grimaldi and Fulde^ there is also a long-range 
contribution, for which this cancellation does not hold, and which is larger than the 
short-range contribution in the long-wave limit. 


The Yafet method is so widely accepted that it may be worthwhile to mention that ours 


is not the hrst paper to point to these problems. For instance. Ref. 


13 


states 


Obviously, [the magnetic moment of state (A, k)] is different for different 
A;-vectors, but it is about ±/iB as long as the spin mixing is small. Yafet neglects 
the /c-dependence of and assumes a constant magnetic moment m or —m 
for all dominant spin-up or all dominant spin-down states, respectively, which is 
certainly problematic in systems with spin “hot spots,” i.e., with regions in the 
Brillouin zone where the spin mixing is very large. 


Also, for semiconductors, Yu et ah— have used spin projection operators to avoid the hrst 
problem in their calculation of the EY spin-relaxation time, but missed the Overhauser 
contribution, as they only include spin-diagonal scattering processes. 

Going beyond spin-relaxation times, there are microscopic approaches to spin-dependent 
carrier dynamics that compute t he reduced density matrix and thus avoid using the pseu¬ 
dospin polarization,=“i£ see Ref. 
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for a review. Numerically solving for the full dynamical 


spin density matrix also yields Dyakonov-Perel and Bir-Aronov-Pikus contributions spin 
relaxation, but is very GPU-time intensive and difficult because the microscopic carrier 
dynamics has to be computed to obtain spin relaxation times that are often orders of mag¬ 
nitude longer than typical scattering times, which have to be resolved in the numerics. 
Numerically calculating the spin density matrix does not yield explicit expressions for spin 
relaxation times, which is one of the goals of our paper. Further, our calculation gives a 
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more transparent description of the ensemble spin dynamics in terms of the torqne matrix 
element than is possible by “brnte-force” calculation of the dynamical spin-density matrix. 

The Elliott-Yafet mechanism has also been applied to the demagnetization dynamics in 
ferromagnets.-'— None of these approaches avoids all of the three problems above. 


III. SPIN-DEPENDENT ELECTRON-PHONON SCATTERING DYNAMICS 

In this section we present a derivation of spin-dependent carrier dynamics due to electron- 
phonon interactions that avoids the three problems listed above. We hrst derive the interac¬ 
tion hamiltonian between electrons and phonons including the long-range contribution of the 
Coulomb potential, and specialize to the long-wavelength limit. With this interaction we de¬ 
rive the equation of motion for the spin density matrix including phonon-assisted correlation 
functions. We identify the torque matrix element that determines the incoherent dynamics 
of the average spin. Finally, we derive the scattering limit of the dynamical equation. 


A. Electron-phonon interaction hamiltonian 

We start our derivation by writing down the interaction potential of a single electron 
in a lattice in terms of lattice site coordinates R„ within the rigid-ion approximation as 
he-L(x, {R„}) = En Wion(x - R„), where 

he-ion(x) = neff(x) - ^ [Vneff(x) X p] ■ S . (8) 

Here, Uefr is the effective electrostatic Coulomb interaction between the electron and the ionic 
core and ^ = h/(4m^c^). We denote single-particle operators acting on the electronic space 
and spin variables by small letters and a hat; for instance, p = —ihV. The second term in ([8]) 
is the spin-orbit interaction. We remind the reader that, due to lack of rotational symmetry, 
this potential does not commute with the electronic angular momentum operator x x p of 
the electrons, and therefore does not conserve the electronic orbital angular momentum. The 
rigid-ion approximation places no fundamental restriction on the following development, but 
the equations of motion become more complicated without it. 

We follow the usual treatment of the electron-phonon coupling by considering small de¬ 
viations R„ = Rn°^ -I- Q„ around the equilibrium conhguration Ri°^ and expand nefr(x — R„) 
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with respect to R„. The resulting electron-phonon hamiltonian in second quantization is 


He-pn = J2 


Q. 


^^e-ion(^ Rn) 


aR, 


R! 


(0) 


T(x) (fx, 


(9) 


where T = (T-]-, denotes a spinor field operator. Note that She-ion/^Rn is an operator 
in spin space. 

While the properties of the electron-phonon interaction can be discussed without specify¬ 
ing a model for the single-particle band structure, we also need the electronic single-particle 
contribution to the hamiltonian for the derivation of the equations of motion. We thus 
consider the model hamiltonian H = He + Re-pn with 

He= /'Tt(x)[^ + ^he-ion(x-R®)]T(x)d3x. (10) 

Jv ^^0 „ 

We follow the standard approach^^ and expand the field operators according to 


where 


T(x) 


1 

TP 


EE 

fi k 


( 11 ) 


‘^pk(x) = (12) 

are Bloch spinor wavefunctions with band and crystal momentum labels (/ik). c^k is the 
corresponding destruction operator, and V is the normalization volume of the crystal. We 
use the convention that a Greek index includes band and pseudospin index, i.e., /i = {b, A). 

Instead of the Bloch wavefunctions, which are orthogonal on the whole crystal volume V, 
we will mainly deal with the lattice periodic n^k’s, which are defined and orthogonal on the 
unit cell G, because we intend to derive matrix elements on the unit cell G. Matrix elements 
of any single-electron operator a that occur in the present paper have the meaning 


(a^kloa^k')!! := E /l“^k(x)l, [atip.k'(x)]^ciV. 


(13) 


By virtue of Eq. ([8]), the electron-phonon interaction hamiltonian ([9]) can be split into 
two terms 


+ Hit"’ ( 14 ) 

the first one of which we call Elliott and the second Overhauser hamiltonian. The Elliott 
hamiltonian itself is spin diagonal 

^^.T = -EEEy/ </5^k(x)Qn • Vnefr(x - R()^^)<y9^/k' (x) Sx, (15) 

n k,k' •' 
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whereas the Overhauser hamiltonian contains the electronic spin operator s 


(x) (Q^ • - R®)) X p] ■ S (x) £x. (16) 

n k,k' 

We next expand the phonon displacement operators according to^^ 

Q„(() = * • (1^) 

q,A 

with Rq ;^ = &LqA + ^q,A; whcrc annihilates (creates) phonons with momentum q in the 
first Brillouin zone and mode index A. We use the abbreviation of = ^yh/{2gV^J^^, 
which has the unit of length, where g is the density. The phonon dispersion is denoted 
by i^q,A- Further, £q a is the polarization vector of the phonon mode with the property 
(£^q,A)* = —^-q,A- We iutroduce the Fourier transformation of the interaction potential in 
the form^ 

= V “=* + G) ( 18 ) 

qelBZ G 

which has a long-range part (G = 0) and short-range part (G 7 ^ 0). 

Now we split the integral over the crystal volume V according to 

y ^ / (x) ^ ^ ^ / (x + d^x, (19) 

use the periodicity M^k(x -|- Rn°^) = M;xk(x) and the relation '-(q+k-k') _ iV5qk-k', 

which neglects Umklapp processes. We thus obtain for the e-pn coupling hamiltonian (I9j) 

He-pn = 5'^k+q,/i'k-®q>^'^/ik+q'^A‘'k ) (20) 

kq /i/i' A 

with the matrix element 

fi'Ajk+q,/i'k ~ ('^Aik+qlW+q,k^M'k)o (21) 

in terms of the e-pn interaction operator 

4+q,k = + G) [£q,A ■ (q + G)] { ^ + (-iO[(q + G)^x (p + hk)] ■ s }. 

G Elliott Overhauser 

( 22 ) 

The “ 1 ” in the curly braces in the above expression comes from the Elliott term (lT5l) and 
the rest from the Overhauser term (IT^ . The q and k vectors are restricted to the first 
Brillouin zone, and we have already neglected Umklapp processes, which contribute if k-f q 





lies outside the first BZ. The long-range part of this expression results from the G = 0 
contribution to the sum, whereas the sum over the G 7 ^ 0 dehnes the short-range part of the 
matrix elements. Grimaldi and Fulde^ demonstrate that the long-range part of the matrix 
element is most important in the long-wavelength limit, so that we explicitly isolate the 
long-range part in the following. 

We will use the e-pn interaction matrix element in the long-wavelength limit as is custom¬ 
ary in semiconductor spintronics.— Grimaldi and Fulde^ also demonstrate how the long- 
range contribution is screened so that neff(q) has a well dehned q —?■ 0 (long-wavelength) 
limit. In this limit, the long-range interaction operator is 

^ 0) (q■ ^q,A){X (p + hk)] -s]. (23) 

Overhauser 

Note that for small q the electrons couple exclusively to longitudinal phonons, where the 
polarization vector points in the same direction as the wave-vector q. 

B. Equation of motion for spin density-matrix 

In this section we derive the equation of motion for the reduced spin density matrix 

P'k' = (cJ,kC/.'k) (24) 

which determines all single-particle properties of the electronic ensemble. 

We will for dehniteness also include coherent terms and eventually take the scattering 
limit. Therefore we need the single-particle hamiltonian ([TU]) in diagonal form 

He = ^ ^ C/.ikc|^k'^Aik- (25) 

Here, we use k ■ p theory,— >2^ i.e., the u^k’s and energy dispersions e^k are determined as 
the solution of the eigenvalue problem of the self-adjoint k ■ p operator heff(k) 

heff(k)'U^k *^/ikn^k (26) 

for lattice-periodic on G. This guarantees that, for each k, the eigenvalue spectrum 
e^k is discrete and the n^k form a complete set, i.e., 

\u^'k){u^k\ = 1. (27) 
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Further, the u^k’s for each k are orthogonal with regard to the scalar product flT^ . In 
semiconductors, the k-p hamiltonian operator can be approximated on a subset of bands by 
an effective hamiltonian, such as the Luttinger hamiltonian.-'^ The number of bands included 
in (123]) determines the dimension of the reduced density matrix pk- 

The equation of motion for the electronic spin-density matrix p|^^ due to (l2Up and (1251) 


is 


with 


Pk -^ifok-VkjPk + ^Pk 


dt 


e-pn 


dt 


Pf 


e-pn 


{ ('*^vk+ql'^k+q,k^Aik)r2(-Bq,ACj^k+q'^Fk) 

A vq 

~ ('W^'k|'l^k,k+q’^!^k+q)o(-B-q,Ac|^k'^i^k+q) }• 


(28) 


(29) 


Equation (|28l) is the full equation of motion that contains the coherent contribution as well as 
the e-pn interaction contribution (129]) . in which the phonon-assisted electronic density matrix 
{^q,Ac|)i^_,_qC^'k) appears.— >2^ The Elliott-Yafet mechanism arises from the e-pn interaction 
contribution in (128]) . To study the EY mechanism in isolation, we will in the following 
neglect the coherent contributions to the equation of motion of the spin density-matrix. 
When we specialize the results to the case of a pair of Kramers degenerate bands in Sec. IIVI 
below, the coherent contributions vanish exactly. 


C. Dynamics of the average spin due to incoherent e-pn interaction 


Before we make approximations to the phonon-assisted density matrix, we compute the 
dynamics of the average spin 




By combining (129]) and 

e-pn n 

A ufifi' k,q 

',F) \ /R J 


T^(x)s„T(x) £x\ = ^ y^(u^k|gaVk)o Pk^'- 

I / 

k 

we obtain 


(30) 


tLs 


(31) 


^ ](K,ik)n(llq,AcJ,,^_|,qC^'k) (^/I'kll'kk+q^kk+q) Q(ll-q,AcJ,,^C^k+q) } - 

We can employ the completeness property (]27]) of the u^k’s to show that this has the form 


m 


Sn 


e-pn 


(^q,A''l'k+q''kk){%'k+ql (tk+q.k)q^kk)£i, 


(32) 


kq 
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with the torque vector operator due to the e-pn interaction 


iw ._lr§y(A) 1 

‘'k+q.k ■- -k ‘^k+q.kj' 


( 33 ) 


We will denote the matrix element of this operator by 


^Ai'k+q,/ik (W'k+q|fk+q,kWk) 


(A) 


(34) 


This matrix element completely determines the spin change that occurs in a transition 
(/ik) (/i'k + q). This information is not explicit the matrix element ([21]). Using (|22|), one 

finds 


tSq,k = + G) [£q,A ■ (q + G)] { ^ + (-iO[(q+ G) x (hk + p)] x s }. 

G Elliott Overhauser 

(35) 

The expression fl32|) together with (l3^ is already an important result, as it shows thatonly 
the Overhauser contribution to the e-pn interaction gives rise to a non-vanishing torque (l35|) . 
The Elliott contribution is spin diagonal and therefore its contribution to the torque vanishes. 
We stress that this result is general enough to include materials with Kramers degeneracy 
(and “spin hot spots”) as well as non-spin degenerate systems with avoided level crossings 
due to spin-orbit coupling. In all cases, the Elliott contribution to the torque matrix element 
does not change the average spin. This conclusion does not depend on the single-particle 
basis in the sense that one could apply a k-dependent unitary transformation to the 
M^k’s and the c^k’s in (132|) and (l35|) . This argument shows that it is not the spin mixing in 
the wave functions that determines the Elliott torque, and—in this restricted sense—there 
is no Elliott contribution to spin dynamics in any basis. 

The torque operator which is due to the Overhauser contribution only, has the 

long-wavelength limit 

tk+q,k = -^^^q°A^eff(q ^ 0) (q ■ £q,A) [q X {Hk + p)] X S. (36) 

Finally, we note also that (j3^ and the definition of the torque matrix element is a rather gen¬ 
eral result for the spin change due to any incoherent electron-boson scattering, for instance, 
with photons or magnons. 
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D. Dynamics of the average spin in the scattering limit 


To compare the results of the present approach with the conventional EY analysis, we 
specialize the dynamical equation (152]) for S'q, for the case of incoherent scattering. To this 
end, we need to evaluate the phonon-assisted density matrix (&q,Ac|jk+q^/i'k) • We do this by 
truncating its equation of motion at the scattering level to obtain 


(^q,Ac|jk+q'^M'k) 


^ ^^tk+q + ^q,A ifl'y) (^q,AcJjk+qC/i'k) 

^ - Pk''Vk+q(l + ^q,A) “ [^r/. 

tt' 



(37) 


where Yq^A = &(^q,A) are phonon occupation numbers, and performing a Markov approxi¬ 


mation, see, e.g.. Refs. 


24 and 


25 


The equation of motion for the phonon-assisted density matrix and, consequently, for 
the reduced density-matrix is determined by the full matrix element including both the 
Elliott and Overhauser contributions, and one finds for the ensemble averaged spin 


dt 


^ S' ? + ift7 

X E - pTipJkd + A'.p) - I'5v - 


(38) 


We will assume that the phonon system is not changed appreciably by the electronic dy¬ 
namics so that the iVq^A = b{hwq^x) are the equilibrium phonon occupations given by the 
Bose function h{huci,x)- At this level of generality, we still have to deal with the reduced 
density matrix on the right-hand side. 


IV. RELAXATION TIME FOR KRAMERS DEGENERATE BANDS 

We now specialize to the case of a pair of Kramers degenerate bands with a quasi¬ 
equilibrium spin-expectation value to obtain a relation for the spin relaxation, which gener¬ 
alizes Yafet’s treatment. 
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A. Eigenstates for Kramers degenerate bands 


Kramers degenerate states including spin-orbit coupling can be assumed to have the 
form2^^— 


Mkir(x) = ak(x)| t) + &k(x)| I), (39) 

Mk^(x) = ahk(x)| ;) - &lk(x)| t)- (40) 


Part of this assumption is that these states are labeled by pseudo spin indices and that, 
generally, ^ |6k(x) P d^x is much smaller than ^ |ak(x) p d^x. Because of the degeneracy, 

there is an ambiguity to the definition of these states, as any superposition will also be an 


eigenstate. In accordance with Refs. 


26 


28 


we choose them to fulfill 


(%,k|s^Mj^,k) = 0. (41) 

This is equivalent to the statement that the chosen states M||.k, M.(j,k diagonalize in the 
degenerate subspace. The choice is important, if one wants to attach some importance to 
the magnitude of ^ \b]^{x.)\‘^ d^x as spin mixing parameter. In ab-initio calculations one 

often uses a small external magnetic field to enforce a quantization direction so that fjdTl) is 
fulfilled. 

The condition (HTl) does not imply that the electrons are in pure spin-up or spin-down 
states characterized by the same two-dimensional spinors at each k. In fact, we have 


"^-fl-jk 7 ^ I T) and 7 ^ | J'); (42) 

as well as |(MAk|'S 2 MAk)| < b,/2 for A =11, 1|. For more details, see appendix R1 


B. Quasi-equilibrium with spin polarization 

In order to derive a characteristic rate for the relaxation of a small excess spin polar¬ 
ization 5Sz at a temperature T, we first need to determine the reduced density matrix for 
this case. We assume the system to be in a quasi-equilibrium with a given spin-expectation 
value 5Sz- Thus we must determine the quasi-equilibrium distribution for the eigenenergies 
of the generalized grand-canonical k • p hamiltonian 

l(k) = heff(k) - Ai - CzSz . (43) 
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Here, the spin accumulation, is the Lagrange parameter needed to obtain the hnite spin 
expectation value. Since heff(k) includes spin-orbit coupling it does not commute with Sz, 
so that the eigenstates and eigenenergies of 

i(k)tc^k = (44) 

are, in general, not identical to the that result as eigenfunctions of /lefr(k), cf. Eq. fl26|) . 

We assume that Cz is small and determine perturbatively in the degenerate subspace 
of the M^k’s. Since we chose the n^k’s that diagonalize Sz in the degenerate subspace, cf. (HIT) , 
this gives 

E^k — Cz I • (45) 

If we also assume that other bands are sufficiently far away, the eigenvectors are w^k — w^k- 
If there are other bands close to the bands for which the spin relaxation is computed, one 
needs to include an eigenvector correction for tc^k- 

The quasi-equilibrium reduced density matrix in the eigenbasis of the grand canonical 
hamiltonian is diagonal and depends only on the grand-canonical energies E^kj he., 

pT' = . ( 46 ) 

where /(e) = [exp(/9e) -|- 1]“^ is the Fermi function. Since (z is small, we keep only the hrst 
order result in (z and with fITS)) we hnd for the quasi-equilibrium distribution 

[/(ek - h) + ^kCz{ui,k\szUf,k )], (47) 

where we have dehned Ak = — For low temperatures, this function approaches a 

5-function concentrated at Ck —/i. From the reduced density matrix (H71) we hnd the relation 
(z = SSz/M between 5Sz and (^z with the normalization 

AT = y^{n^k|szM/.k)^Ak . (48) 

/ik 


C. EY relaxation time 


To obtain the close-to-equilibrium spin relaxation rate dehned by 


d 


SSz = 


SSz 


(49) 


dt tsr 

we insert (lT7j) in fl38)) and linearize. Then, as in Yafet’s original treatment and its extension 
by Fahnle and coworkers,— (z drops out and we hnd 
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( 60 ) 


(^t'k+q,Aik)^5'^k,;j'k+q('*^Mk+q|-S^M/xk+q)Ak+q 

kq fifi' A 


X 


-^q,A “1“ f (*^/^k /^) 


+ 


1 + ^q,A — /(C/^k — A'-) 


■ ^^i'k+q *^/^k “1“ ^^^q,A “1“ (^fi'k+q ^^q,A “1“ ■ 

This expression for the spin relaxation, or Ti, time is a more physically transparent result 
compared to Yafet’s because it relates the spin dynamics directly to a torque matrix element. 
It improves on the original Yafet result by including the correct bookkeeping for spin instead 
of accounting for pseudo-spin flips. Instead of P; i-e-, the squared modulus of the e-pn 

matrix element, as it occurs in Yafet’s result, here matrix elements of the torque, the electron- 
phonon interaction and the spin appear. As shown in Sec. IlllCI the torque matrix element 
only has an Overhauser contribution, whereas in the g matrix element both Elliott and 
Overhauser terms contribute. Regarding the g matrix element in the long wavelength limit, 
Yafet^ and numerical evaluations^^ of his spin-relaxation result find a cancellation between 
Elliott and Overhauser contributions in the short range contribution, whereas Grimaldi and 
Fulde^ find for the long-wavelength limit a dominating long-range Elliott contribution to g. 

Our result fl50|) is not limited to small spin mixing and is qualitatively different from the 
original Yafet formula, as there is no pure Elliott contribution. It will make it easier to com¬ 
pute accurate EY spin relaxation times numerically and to compare different contributions 
to spin dynamics by computing the relevant combination of torque and g matrix elements. 


V. CONCLUSION 

In conclusion, we have presented an investigation of the carrier-spin dynamics due to 
incoherent electron-phonon scattering in the presence of spin-orbit coupling, as originally 
considered by Overhauser, Elliott, and Yafet. We examined the dynamical equation for 
the reduced density matrix in terms of the phonon assisted density matrix including the 
contributions of the spin-diagonal (electrostatic) electron-phonon interaction as well as the 
phonon-modulated spin-orbit coupling. We showed that the central quantity to account 
for the spin change is a torque matrix element, which gives a physically more appealing 
picture of spin dynamics due to the EY mechanism. Our results show that the electron- 
phonon “Elliott contribution,” i.e., the spin mixing in the Bloch states, has no impact on the 
torque matrix element due to the electron-phonon interaction. This important property of 
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spin-dependent carrier dynamics gets lost if one works only with pseudo-spin distributions 
because these do not characterize the spin completely. Only the explicitly spin-dependent 
phonon-induced modulation of the spin-orbit coupling as introduced by Overhauser gives rise 
to a non-vanishing torque due to an incoherent scattering mechanism. Finally, we presented 
an explicit expression for the spin relaxation, or Ti, time close to equilibrium, which is valid 
for arbitrary spin mixing. 
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Appendix A: Spin expectation values for degenerate bands 


We show here in detail that, for the case of degenerate bands, the electron spin operator 
can be diagonalized in the space of the degenerate bands, even if spin-mixing due to spin- 
orbit coupling is present. This is not a new result, but it is important to keep this in mind 
when dealing with Kramers degenerate wave functions. 

We assume we have arbitrary tick’s of the general form 

Mik(x) = ak(x)| t) + &k(x)| I), (Al) 

M2k(x) = ahk(x)| ;) - &hk(x)| t), (A2) 


which are Kramers degenerate. In this subspace, for hxed k the matrix of the spin operator 
is 


h I dii 


Sz = 



(A3) 


with the real diagonal elements 


dk • ("^kl 1^2:1 "^kl) 



(A4) 


and complex off-diagonal elements 

Ck • ('^k2 1^2 I "^kl) 



[ok (r) 6_k (r) + a-k (r) &k (r)! • 


(A5) 
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If Ck 7 ^ 0, diagonalization leads to the eigenstates maic with A ='||, 


/ \ Ck*| dk “1“ d_k 

%k(x) = J I t) - 

dk “1“ d_k 


1 + 


4 |ck| 


(dk + d_k) 


%k(x) = ^1;) + 


2 L 


1 + 


4 |ck| 


(dk + d_ 


U), 

t), 


where L normalizes the maic’s. The corresponding spin eigenvalnes are 

h 


X 


^(dk - d_k) ± Y ^(dk + d_k)2 + |ckp 


(A6) 

(AT) 


(A8) 


In the presence of inversion symmetry, as in alnmininm and silicon, these expressions simplify 
fnrther. Then dk = d_k and the eigenvalnes are ±(h/2)-^d^ + |ckP for A =ft', IJ-, respectively. 
This resnlt shows that psendo spin labels are jnstihed as these states generally have an 
expectation valne of Sz of modnlns smaller than (or almost eqnal to) ±h/2. At the same 
time, the psendo-spin states are not jnst rotated pnre spin states, which wonld lead to a 
spin projection of exactly ±h/2 along the qnantization axis. Fnrther, these states may have 
nonvanishing and Sy expectation valnes. 
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